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Henry  F.  Schaefer  III 


Fundamental  Studies  of  Carbon,  NH,  and  Oxygen  Rings  and  Other  High 
Energy  Density  Molecular  Systems 


The  development  of  efficient  and  safe  conventional  (i.e.,  nonnuclear)  propellants 
and/or  fuels  is  a  goal  of  obvious  technological  significance.  A  desirable  quality  of 
such  a  propellant  is  clearly  a  high  ratio  of  energy  release  to  mass.  The  present 
hypothesis  rests  on  a  simple,  but  previously  unrecognized,  analogy  between  oxygen 
and  sulfur.  Preliminary  studies  showed  that  the  oxygen  ring  systems  are  sufficiently 
promising  to  warrant  the  detailed,  high-level  theoretical  research  reported  here.  -  - 

Our  idea  begins  with  the  observation  that  elemental  sulfur  exists  as  sulfur  rings, 
Sn.  The  stability  of  these  sulfur  rings  is  attested  by  the  fact  that  the  single  most  impor¬ 
tant  commercial  source  of  sulfur  is  elemental  sulfur  in  the  caprock  salt  domes  in  the 
United  States  and  Mexico.  The  most  common  allotrope  of  sulfur  is  the  yellow, 
orthorhombic  a  form,  to  which  all  other  modifications  eventually  revert  at  room  tem¬ 
perature.  As  long  ago  as  1935  the  crown  structure  of  cyclo-S8  was  established  by  x-ray 
crystallography.  The  D4d  point  group  experimental  structure  of  S8  is  widely  known. 

The  essence  of  our  proposal  is  to  make  an  analogy  between  sulfur  rings  and  oxy¬ 
gen  rings.  Given  the  remarkable  stability  of  sulfur  rings,  should  it  not  be  possible  to 
prepare  oxygen  rings?  Oxygen  lies  directly  above  sulfur  in  the  Periodic  Table,  and  the 
analogy  is  an  appealing  one.  Before  pressing  on  further  with  this  analogy,  let  us  say  a 
bit  more  about  cyclic  forms  of  sulfur  other  than  S8. 

Although  S8  dominates  elemental  sulfur,  more  than  a  dozen  additional  sulfur  rings 
Sn  have  been  synthesized  in  the  past  20  years.  Long  before  this,  in  1891,  a  rhom- 
bohedral  form,  e -sulfur,  was  prepared  by  Engel.  By  1961  € -sulfur  was  clearly  esta¬ 
blished  to  be  cyclo- S6.  Note  that  for  both  S6  and  S8  all  sulfur  atoms  are  equivalent.  S7 
is  alos  known  to  exist  as  a  Cs  ring  structure,  with  four  distinct  S-S  distances,  ranging 
from  1.99  to  2.18  X.  Sn.  and  S13  have  been  synthesized  and  appear  from  spectroscopic 
studies  to  be  cyclic,  although  no  crystal  structures  have  been  reported. 

S10  belongs  to  the  rarely  observed  point  group  D2  and  has  experimental  bond  dis¬ 
tances  in  the  relatively  narrow  range  2.03  -  2.08  X.  The  12  membered  sulfur  ring  was 
predicted  in  1949  in  a  classic  paper  by  Pauling  to  be  unstable.  Nevertheless,  S12  was 
synthesized  in  1966  and  turns  out  to  be  second  only  to  S8  in  stability.  S12  is  monocy¬ 
clic,  with  S-S  bond  distances  in  the  very  narrow  ra-^e  2.C5  -  7  06  X.  S28  is  also  mono- 
cyclic,  as  is  S^  and  experimental  cyrstal  structures  exist  for  both  molecules. 
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The  reason  for  this  relatively  detailed  description  of  the  cyclic  Sn  isomers  is  that 
there  should  be  a  one-to-one  correspondence  with  cyclic  On  isomers,  a  principal  topic 
of  this  research.  However,  unlike  the  Sn  sulfur  rings,  the  analogous  On  rings  will  be 
high  energy  density  molecular  systems. 

Before  going  on  further,  let  us  attempt  to  make  some  plausible  guesses  as  to  the 
energy  content  of  the  metastable  On  rings.  First  one  notes  that  the  valence-isoelectronic 
cyclic  S8  contains  no  energy  in  this  sense.  Specifically,  crystalline  S8  lies  15.4  ±0.1 
kcal  per  mol  of  sulfur  atoms  below  four  gas  phase  diatomic  S2  molecules.  Can  one 
guess  thsi  result  using  the  dissociation  energy  of  S2  (D0  =  4.37  eV  =  100.8  kcal/mol) 
and  some  reasonably  standard  S-S  single-bond  dissociation  energy,  say  54  kcal/mol? 
The  answer  to  this  question  is  a  qualified  "yes".  Using  this  simple  model,  one  predicts 
S8  to  lie  (54  -  50.4)  =  3.6  kcal/mol  below  four  S2  molecules  on  a  per  atom  basis.  By 
increasing  the  S-S  bond  energy  from  the  standard  54  to  66  kcal,  the  known  experimen¬ 
tal  energy  difference  for 

Sg(crystalline)— »S2  1 

is  precisely  reproduced.  This  adjustment  reflects  the  fact  that  the  S-S  bonds  in  Sg  are 
stronger  than  those  in  organosulfur  compounds  such  as  CH3SSCH3, 

The  comparable  oxygen  thermodynamic  data  make  it  immediately  obvious  why 
oxygen  rings  should  be  high  energy  density  materials.  First  the  standard  0-0  bond  dis¬ 
sociation  strength  is  35  kcal/mol,  much  weaker  than  the  54  kcal  for  S-S  bonds. 
Second,  the  dissociation  energy  of  diatomic  02  is  D0  =  5.12  eV  =  118.0  kcal/mol, 
much  stronger  than  the  100.8  kcal  for  S2.  Thus,  the  estimate  for  the  energy  stored  in 
the  "generic"  oxygen  ring  is  (35-59_  =  -24  kcal  per  mole  of  oxygen  atoms.  That  is,  the 
dissociation  process 


On  (n  even) 


f  0,(6) 


is  estimated  to  be  exothermic  by  24  kcal  per  mole  of  oxygen  atoms.  Alternately  48 
kcal  of  energy  is  released  per  mole  of  02  molecules.  A  larger  energy  release  might  be 
expected  for  the  smaller  rings,  specifically  O4  and  06,  which  presumably  have  smaller 
average  0-0  bond  energies  due  to  ring  strain. 

A  similar  estimate  is  readily  made  for  the  isoelectronic  cyclic  system  (NH)n, 
which  will  also  be  metastable  and  dissociate  (directly  or  indirectly)  via 


(NH)n  (n  even)  ->  (g)  +  ~  H2  (g). 


This  proposal,  of  course,  takes  advantage  of  the  tremendous  strength  (D0  =  9.76  eV  = 

255.1  kcal/mol)  of  the  NsN  triple  bond.  Using  in  addition  D0(H2)  =  4.48  eV  =  103.3  oa/ 
kcal  and  the  average  N-N  single  bond  strength  (39  kcal)  and  N-H  bond  strength  (93  .  ty  Codft3 


...ii.  and/or 
Special 


□a 
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39  +  93  -  112.5  -  51.6  =  -32.1  kcal 

per  mole  of  N  (or  H)  atoms.  Alternately,  reaction  (3)  is  suggested  to  be  exothermic  by 
64.2  kcal/mol  of  N2  molecules  releases.  It  is  apparent  that  within  the  limits  of  the  sim¬ 
ple  estimates  made  here,  systems  like  (NH)8  may  be  even  more  effective  than  the 
above-discussed  oxygen  rings  as  high  energy  density  molecular  systems.  It  is  ulti¬ 
mately  intended  to  theoretically  investigate  both  On  and  (NH)n.  The  order  of  the  dis¬ 
cussion  presented  above  was  mandated  by  the  analogy  (with  the  sulfur  rings)  by  which 
this  idea  came  to  us. 

A  traditional  measure  of  the  effectiveness  of  a  propellant  is  the  specific  impusle 
Isp-  TsP  given  in  seconds  from  the  relationship 

T  ,  s  AH(kcal/mol) 

Isp  =  (constant)— — - — — - — — . 

v  Molecular  weight(g/mol) 

The  traditional  standard  of  comparison  is  the  exothermic  reaction 

H2  +  '/202  ->  H20  4 

for  which  Isp  is  on  the  order  of  400  s.  Substitution  of  unstrained  (e.g.,  Og)  oxygen 
rings  for  02  will  increase  the  specific  impulse  by  -100  s.  Use  of  the  smaller  oxygen 
rings  04,  05,  or  06  would  presumably  result  in  a  significantly  larger  Isp. 

It  should  be  emphasized  of  course  that  high  energy  density  is  a  necessary  but  cer¬ 
tainly  not  sufficient  condition  for  the  development  of  an  effective  propellant.  The 
material  also  must  have  a  significant  range  of  metastability  -  it  should  not  explode 
when  a  match  is  unintentionally  lit  100  m  away.  The  great  stability  of  elemental  sulfur, 
in  the  form  of  sulfur  rings,  gives  some  encouragement  that  the  band  of  metastability 
for  oxygen  and  NH  rings  might  be  rather  broad.  Precisely  this  point,  of  course,  may  be 
investigated  in  the  theoretical  studies  of  the  type  reported  here. 

Due  to  the  weakness  of  the  0-0  single  bonds,  as  in  hydrogen  peroxide,  it  is 
hypothesized  that  oxygen  rings  are  potential  high  energy  density  materials.  A  particu¬ 
larly  attractive  candidate  is  the  04  molecule,  for  which  ring  strain  is  expected  to  pro¬ 
vide  further  destabilization  relative  to  two  separated  02  molecules.  To  pursue  these 
qualitative  suggestions,  ab  initio  molecular  quantum  mechanics  has  been  employed. 
Both  self-consistent  field  (SCF)  and  configuration  interaction  including  single  and  dou¬ 
ble  excitations  (CISD)  methods  have  been  employed  in  conjuction  with  double  zeta 
plus  polarization  basis  sets.  At  the  highest  level  of  theory  the  nonplanar  (DM  point 
group,  O-O-O-O  torsional  angle  25  °)  equilibrium  structure  is  predicted  to  lie  2.9  kcal 
below  the  planar  D4h  structure,  which  is  a  transition  state.  The  infrared  spectrum  is 
predicted  at  the  DZ  +  P  CISD  level,  as  well  as  lower  levels  of  theory.  The  04 
minimum  is  predicted  to  lie  -  100  kcal/mol  above  the  asymptotic  limit  of  two  02 
molecules. 

Self-consistent  field  (SCF)  methods  with  minimum  (STO-3G),  double  £  (DZ),  and 
double  £  plus  polarization  (DZP)  basis  sets  predict  the  06  ring  to  assume  chair,  twist, 
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Self-consistent  field  (SCF)  methods  with  minimum  (STO-3G),  double  C,  (DZ),  and 
double  C  plus  polarization  (DZP)  basis  sets  predict  the  0$  ring  to  assume  chair,  twist, 
and  boat  conformations  analogous  to  similar  forms  for  cyclohexane.  All  predicted 
vibrational  frequencies  for  the  chair  and  twist  forms  are  real.  Six  symmetrically 
equivalent  oxygen  atoms  are  predicted  to  comprise  the  lowest  energy  chair  form,  with 
0-0  bond  distances  of  1.364  X  and  bond  angles  of  104.7°  at,  the  DZP  SCF  level  of 
theory.  The  boat  form  is  not  found  to  be  an  energy  minimum  but  rather  exhibits  one 
imaginary  vibrational  frequency  which  when  followed  tends  toward  assumption  of  the 
twist  form.  Energy  differences  at  the  DZP  SCF  level  are  computed  to  be  15.9  kcal 
between  the  chair  and  twist  forms  and  17.5  kcal  between  the  chair  and  boat.  We  inter¬ 
pret  these  results  by  analogy  with  cyclohexane  and  assign  the  larger  energetic 
discrepancies  to  shorter  bond  distances  and  inherently  greater  eclipsing  effects  for  adja¬ 
cent  lone  electron  pairs  than  those  attributed  to  bonding  electron  pairs.  Homodesmotic 
and  hyperhomodesmotic  reactions  devised  to  predict  the  decomposition  exothermicity 
of  the  ring  give  rather  different  results,  namely,  -130  (homodesmotic)  and  -75 
(hyperhomodesmotic)  for  the  heat  for  formation  of  06.  Higher  level  theoretical  studies 
will  be  required  to  resolve  this  energetic  discrepancy. 


